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Abstract 



The Heisenberg chain with a weak link is studied, as a simple example of a 
quantum ring with a constriction or defect. The Heisenberg chain is equiva- 
lent to a spinless electron gas under a Jordan- Wigner transformation. Using 
density matrix renormalization group and quantum Monte Carlo methods we 
calculate the spin/charge stiffness of the model, which determines the strength 
of the 'persistent currents'. The stiffness is found to scale to zero in the weak 
link case, in agreement with renormalization group arguments of Eggert and 
Affleck, and Kane and Fisher. 

PACS Indices: 73.23.Ra, 75.10.Jm, 68.66.La 
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I. INTRODUCTION 



Technological advances in recent years have allowed the fabrication of electrical and 
even mechanical devices on the nanometer scale, where individual atoms or electrons can be 
manipulated. The physics of these devices poses a plethora of fundamental questions through 
a rich variety of novel quantum effects IJ. This has led to an upsurge of theoretical interest 
in the physics of 'quantum wires', 'quantum dots' and more general physics at the mesoscopic 
or nanometer scale |2| . The effects of electron-electron interactions are typically enhanced in 
systems of reduced dimensionality, leading to nonperturbative effects, such as the breakdown 
of Fermi liquid behavior in one- dimensional metals, and single-electron charging effects in 
quasi- zero-dimensional systems (quantum dots). An important milestone in the nascent 
field of nanomechanics was the experimental discovery that not only the electrical but also 
the mechanical properties of metallic structures on the nanometer scale exhibit apparently 
universal nonmonotonic quantum corrections which could be explained theoretically 

within the framework of a Jellium model |J. 

In the electrical domain, paradigm systems to investigate mesoscopic behavior have long 
been small ring-shaped or multiply connected devices, where the application of a magnetic 
flux piercing the device leads to persistent currents. However, the theoretical prediction of 
these persistent currents, first in superconducting and then in normal conducting materials, 
has always predated experimental investigation, which has only become feasible in the last 
decade 0. Very recently, with the discovery of a tunable Kondo effect in quantum dots, 
the persistent currents of multiply connected systems with magnetic quantum dots - Kondo 
rings for short - have aroused considerable interest 10^1 • 

In this paper we explore a very simple system, which may serve to model a metallic 
quantum wire ring with a weak junction, or constriction. It consists of the standard spin- 
1/2 Heisenberg antiferromagnetic spin chain, which by a Jordan- Wigner transformation 
is equivalent to a spinless electron gas in one dimension, where the exchange coupling is 
weakened at a single link. We use Density Matrix Renormalization Group (DMRG) |pT0|JTT 



and Quantum Monte Carlo (QMC) methods [E| to obtain numerical results on chains of 
up to 256 sites, and perform finite-size scaling extrapolations to the bulk limit. We study 
the spin stiffness, which under the Jordan- Wigner transformation is equivalent to the charge 
stiffness of the electron gas and is related to the persistent current, as a function of the weak 
link coupling. We also study the spin correlations across the weak link. 

Eggert and Affleck [|13] have previously studied the Heisenberg chain with an isolated 
impurity using exact diagonalization and conformal field theory techniques. They find that 
in renormalization group language a single weak link across the ends of an open chain corre- 
sponds to an irrelevant operator, and therefore the open chain is a stable fixed point under 
such a perturbation. Thus they predict that in the bulk limit a chain with a weak link will 
behave like an open chain. These findings for a concrete model system are in agreement with 
the general predictions of Kane and Fisher |14| for the general one-dimensional interacting 
electron gas, i.e. the Luttinger liquid. An integrable version of the Heisenberg spin chain 
with defects has also been studied |T5|-p^. The case of a single defect corresponds to a 
weak link. However, as opposed to the case discussed here where only one bond is modified, 
integrability requires a modification of two adjacent bonds and an additional three-spin 
coupling. Although the defect of this integrable chain is completely transparent to particle 
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scattering, the persistent current is renormalized by the defect strength. 

Let us briefly review [18] how the persistent current arises, for the simple case of free 
electrons. Start from the real-space continuum Hamiltonian 

H = -^-E [ ' dx^ a {x)dlU*) (1) 



where ip a (x) is an electron field with spin index a = ±1, and L is the circumference of the 



ring. Now thread the ring with a magnetic flux $, producing an Aharonov-Bohm effect ||T9 
The quantum phase 



L e 



He Jo A ^ X = Tc^ (2) 
can be encoded via a gauge transformation in the twisted boundary conditions 

ML) = e^a(O), (3) 

where 

= 2vr^ (4) 

and $o = hc/e is the elementary flux quantum. 

When an Aharonov-Bohm field is applied, the Hamiltonian acquires the usual interaction 
term (we set % = c = 1 henceforth) 

H int = - J dxA^JHx). (5) 

Thus for a constant field A\ = $/L the corresponding "persistent current" is given by the 
Feynman-Hellman theorem: 



which can be expanded 



m) = -D c j + 0({jf) (7) 



where D c is the "charge stiffness" . If we assume that /($) is purely linear in $ (as can be 
proved for the pure Heisenberg chain |20fl), then the charge stiffness and hence the persistent 
current can be estimated from the difference in energy between the system with anti-periodic 
boundary conditions ( Eq = E (<& = $o/2)) and periodic boundary conditions ((Eq = 
£ ($ = 0)) 



D c = 8 -±iE -E+] (8) 



The corresponding quantity in the Heisenberg chain is the spin stiffness. 

In section II of the paper we briefly summarize the DMRG and QMC methods used to 
calculate this quantity. In section III we present our results, and in section IV our conclusions 
are summarized. 
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II. METHOD 



We study the spin-1/2 Heisenberg quantum spin chain with a single weak coupling J' < J 
between two adjacent spins located between sites i = N and i — 1. The Hamiltonian is 

JV-l 

H = J ^ • S i+1 + J'S N ■ Si. (9) 

8=1 

There are a total of N sites in the ring. We study J'/ J in the range [0, 1], hence we have 
either open or periodic boundary conditions (OBC or PBC) at the extremities of this range. 
We also consider anti-periodic boundary conditions in the same range of J' / J, to obtain the 
charge stiffness D c according to (pi). The model in these limits is exactly solvable by Bethe 



ansatz [21J , and hence is often used as a testing ground for various DMRG methods |10| , |22 
The quantities that have been calculated using DMRG include the ground state energy, the 
singlet and triplet gaps, and correlation functions. 



The "infinite-lattice" DMRG method ITUj] is used here, applied with periodic boundary 



conditions. The lattice is split into two blocks and two sites as shown in Fig. [T]. The weak 
link is placed between block 1 and site 2. At any time the superblock consists of a system 
block and an environment block, plus two extra sites. The presence of the weak link destroys 
the translational invariance normally exploited in usual DMRG schemes, hence we cannot 
simply make a copy of the density matrix in one block and transfer it to the other block. 
Therefore in a single DMRG iteration, two density matrices are constructed (one for each 
block), and the basis set for each block originates from its corresponding density matrix. 
Each block increases in size by a single site in a single DMRG iteration. We calculate results 
for lattice sizes N = 4 to 64, in steps of two. The quantities calculated here are the ground 
state energies, and the correlation of the spins across the weak link. The total number of 
density matrix eigenstates retained in a block was m = 350 in the basis truncation procedure. 

We also carried out quantum Monte Carlo simulations using the stochastic series expan- 
sion (SSE) method jl2| . In this case, the spin stiffness can be directly calculated as the 



second derivative of the energy with respect to the phase, which is given in terms of the 
fluctuation of the winding number in the simulations |f23| . 



III. RESULTS 

To estimate the accuracy of the DMRG, we perform convergence tests with m, the 
number of basis states retained in a block. Table | shows sample results for J' / J = 0.5 and 
N = 64 with both periodic and anti-periodic boundary conditions. We see good convergence 
with the number of basis states retained: in particular, for the ground state energy we have 
convergence to 1 part in 10 6 . We obtained independent estimates using SSE techniques 
for the periodic case, which yields E /J = —28.2178(3), agreeing perfectly with the DMRG 
results. The correlations between the spins across the links have a marginally lower accuracy 
because of round-off error, but even here we have an error in the region of 1 part in 10 4 . 

Since the SSE method is a finite-temperature quantum Monte Carlo method we have 
to run the simulations at sufficiently low temperature to converge the quantities of interest 
to their ground state values. In Table [TJ we show the convergence of the energy and the 
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charge stiffness for a 256-site chain. All SSE results discussed below were obtained at inverse 
temperatures (3 where the results do not differ, within statistical errors, from results at (3/2. 

Ground state energies were calculated for the periodic and anti-periodic rings, as a func- 
tion of the weak link coupling J' / J. Figs. H and |3] shows DMRG estimates for the quantity 
AEn = (Eo(N; J') — E (N; J' = J))/ J as a function of J'/ J, for several different lattice 
sizes up to N = 64. An extrapolation in 1/N can be performed to extract the bulk limit 
for each value of J' by a simple polynomial fit to the data, giving us the extrapolated curve 
for the periodic and anti-periodic cases. It can be seen that for the periodic case the bulk 
values are approached from above, while for the anti-periodic case the limit is approached 
from below. 

Putting together both the periodic and anti-periodic results for the energy, we can cal- 
culate the spin stiffness factor, given by 

2N 

p s = -——(E (N; anti-periodic) — E (N; periodic)) (10) 

TT 2 J 

The results are shown in Fig. ^. At couplings other than J' /J = 1 the values trend steadily 
down towards zero as the lattice size N increases. There is a marked difference in behavior 
for the isotropic case J'/ J = 1, as the strong curvature towards zero is not apparent. One 
would naively expect to obtain a bulk limit by a simple linear extrapolation procedure, 



but in fact the work of Woynarovich and Eckle |24j has shown that there are logarithmic 
corrections to the ground state energy, and hence the stiffness. We can see the effects of 
these corrections for the J' /J = 1 there exists an exact result obtained by Hamer, 



Quispel and Batchelor [E2IJ (equation (3.37) of ref. [^IJ with 7 = 0) 

Ps = ~- (11) 

This does not seem in accord with the data in Fig f|, which appear to be approaching 0.27. 
A late, logarithmic downturn must therefore occur at very large lattice sizes. 

Using field theoretical methods, Eggert and Affleck |13| have predicted that a chain with 
J' < J should be similar to an open chain (i.e., J' /J = 0). This implies that the spin 
[charge] stiffness should vanish as the system size N — > oo. Fig. |5| shows SSE results [25 



for the stiffness versus the system size (N = 16, 32, 64, 128 and 256) for several values of 
J' < J. The results are in accordance with a scaling behavior 

p N (x) ~ a(x)N- u (12) 

where x = J'/ J, and the index a ~ 2/3. Fig. |6] demonstrates that the data for large N can 
be well described, in fact, by a simple scaling form 

9 (St 

Pn(x) ~ jf^N~ 2/3 , N - oo. (13) 

It is likely, however, that the true asymptotic correction-to-scaling behavior is again being 
disguised by logarithmic corrections. 

We have also calculated the value of the spin correlation function across the weak link, 
i.e. (SffS*). Fig. [7| shows the behavior versus 1/N 2 for various values of J' /J. It can 
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be seen that the finite-lattice values generally approach a finite value in the bulk limit, as 
one would expect, except for the special case J'/ J = where the link is open. Theoretical 
expectations [^6j are that the correlation function should approach its bulk limit like l/N 2 , 
up to logarithmic corrections. 

The presence of a finite correlation across the weak link is just what one would naively 
expect when the weak-link coupling J' is non-zero. On the other hand, it might appear to 
contradict the previous statement that a chain with a weak link should renormalize to the 
open chain. The point here is that the weak-link correlation is a local quantity, not a bulk 
property. It is only bulk properties such as the spin-stiffness which scale to the value of the 
open chain. 



IV. CONCLUSIONS 

In summary, we have performed a finite-lattice study of the Heisenberg ring with a weak 
link, using both DMRG and QMC calculations on rings of up to 256 sites. The spin or 
'charge' stiffness has been calculated either directly (QMC), or from the energy difference 
between the system with anti-periodic boundary conditions and that with periodic bound- 
aries, assuming a quadratic dependence of the energy on the twist parameter <fi (DMRG). 

The stiffness, and hence the persistent current, is found to scale to zero in the bulk limit 
iV — ► oo, for any J' < J. This agrees with the renormalization group prediction of Eggert 
and Affleck |13] , that the stable fixed point for this system corresponds to an open chain, so 
that the chain with a weak link will behave like an open chain, as regards its bulk properties. 

We have also measured the spin-spin correlation across the weak link. A finite antiferro- 
magnetic correlation remains in the bulk limit, depending on the coupling J' as one would 
expect. The renormalization group argument does not apply to a 'local' quantity such as 
this. 

The asymptotic scaling behavior of these quantities has been disguised to some extent 



by logarithmic finite-size scaling corrections. Eggert and Affleck [13] have circumvented this 
problem by studying a modified model with an extra marginal operator; but we have not 
found this necessary for our present purposes. 

For the future, it would be of interest to see how the results generalize to more compli- 
cated and interesting cases, such as higher spin chains, or real electronic models, such as the 
Hubbard model or its variant, the so-called t — J model. Another interesting extension of 
the present study would be to interpret the weak link and hence the modified bond in our 
model as caused by a mechanical force on a quantum wire. It would be interesting to see 
what conclusions could be drawn from our simple one- dimensional model for such a scenario. 
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TABLES 



TABLE I. DMRG estimates of the ground state energy Eq/J and correlation (S^Sf) for 
N = 64 sites at J'/J = 0.5 as a function of m, the number of states retained per block. A 
SSE estimate of Eq/J = -28.2178(3) for PBC agrees very well with the DMRG data: here our 
final DMRG estimate is -28.21797(1). The correlation between spins across the weak link also 
converges to better than 1 part in 10 4 . The anti-periodic boundary conditions yield similar levels 
of accuracy. 



m 


Periodic 

E /J 


(S z N Sf) 


Anti-periodic 

E /J 


(S z N Sf) 


96 


-28.217652 


-0.060888 


-28.212779 


-0.0560186 


164 


-28.217938 


-0.061194 


-28.213082 


-0.0562742 


234 


-28.217964 


-0.061219 


-28.213123 


-0.0563028 


342 


-28.217970 


-0.061222 


-28.213132 


-0.0563091 



TABLE II. The internal energy and the spin stiffness calculated in SSE simulations for a 
256-site chain with J'/J = 1/4 at different inverse temperatures = J/T. 



J3 E /NJ Ps 

32 -0.44211(2) 0.0000(0) 

64 -0.44237(1) 0.0008(1) 

128 -0.442429(9) 0.0088(3) 

256 -0.442443(5) 0.0202(5) 

512 -0.442443(4) 0.0212(4) 

1024 -0.442453(3) 0.0215(2) 

2048 -0.442448(3) 0.0218(2) 
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FIG. 1. The augmentation process within one DMRG iteration. Augmentation 1 (Augmenta- 
tion 2) gives the new Block 1 (Block 2) in the next DMRG iteration. 
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FIG. 2. AE N = (E (N; J') - E (N; J' = J))/ J as a function of J'/J for the ring with periodic 
boundary conditions. The data is extrapolated using a simple fit to obtain the bulk limit. 
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FIG. 3. As for Fig. ||, but with anti-periodic boundary conditions. 
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FIG. 4. The stiffness factor p s as a function of 1/N, for lattice sizes N = 4 to 64. 
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FIG. 6. Scaling plot of p s /x versus (1 — x) 3 ^ 2 N, where x = J' /J. Also shown is the scaling 
form (O), which agrees with the data for large N. 
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FIG. 7. The correlation (Sf^Sf) across the weak link as a function of 1/N 2 
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